MATRIX P-NORMS ARE NP-HARD TO APPROXIMATE IF 

P ^ l,2,oo. 

JULIEN M. HENDRICKX AND ALEX OLSHEVSKY * 

Abstract. Wc show that for any rational p G [1, 00) except p = 1, 2, unless P = NP, there is no 
polynomial-time algorithm which approximates the matrix p-norm to arbitrary relative precision. Wc 
also show that for any rational p £ [1, 00) including p = 1, 2, unless P = NP, there is no polynomial- 
time algorithm which approximates the 00, p mixed norm to some fixed relative precision. 

1. Introduction. The p-norm of a matrix A is denned as 

\\A\\ p = max ||Ab|| p . 

||x|| p =l 

We consider the problem of computing the matrix p-norm to relative error e, defined 
as follows: given the inputs (i) a matrix A <E R nxn with rational entries (ii) an error 
tolerance e which is a positive rational number, output a rational number r satisfying 

\r-\\A\\ p \<e\\A\\ p 

We will use the standard bit model of computation. When p = 00 or p = 1 the p- 
matrix norm is the largest of the row/column sums, and thus may be easily computed 
exactly. When p = 2, this problem reduces to computing an eigenvalue of A T A and 
thus can be solved in polynomial time in n, log - and the bit-size of the entries of A. 
Our main result suggests that the case of p ^ {1, 2, 00} may be different: 

Theorem 1.1. For any rational p £ [1, oo) except p = 1,2, unless P = NP, there 
is no algorithm which computes the p-norm of a matrix with entries in {—1,0, 1} to 
relative error e with running time polynomial in n, i. 

On the way to our result, we also slightly improve the NP-hardness result for the 
mixed norm ||vl||oo,p = max ||:c||oo<i. II^IIp from [5]. Specifically, we show that for 
every rational p>l, there exists an error tolerance e(p) such that unless P = NP, 
there is no polynomial time algorithm approximating H-AHo^j, with a relative error 
smaller than e{p). 

1.1. Previous work. When p is an integer, computing the matrix norm can be 
recast as solving a polynomial optimization problem. These are known to be hard to 
solve in general [3j; however, because the matrix norm problem has a special structure 
one cannot immediately rule out the possibility of a polynomial-time solution. A few 
hardness results are available in the literature for mixed matrix norms ||y4.|| p ,g = 
max 1 1 j. 1 1 <! HArllq. Rohn has shown in 0] that computing the ||^4||oo,i norm is NP- 
hard. In her thesis, Steinberg [5] proved more generally that computing ||^4|| Pjq is 
NP-hard when 1 < q < p < oo. We refer the reader to [S] for a discussion of 
applications of the mixed matrix norm problems to robust optimization. 
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It is conjectured in [5] that there are only three cases in which mixed norms are 
computable in polynomial time: p = lorg = ooorp = g = 2. Our work makes 
progress on this question by settling the "diagonal" case of p = q; however, the case 
of p < q, as far as the authors are aware, is open. 

1.2. Outline. We begin in Section [2] by providing a proof of the NP-hardness 
of approximating the mixed norm || • ||oo, P within some fixed relative error, for any 
rational p > 1. The proof may be summarized as follows: observe that for any matrix 
M, max|| x || oo=1 ||Mx|| p is always attained at one of the 2" points of { — 1, 1}"; so by 
appropriately choosing M , one can encode the NP-hard problem of maximization over 
the latter set. This argument will prove that computing the || • ||oo, P norm is NP-hard. 

Next, in Section [3] we exhibit a class of matrices A such that max|| a ,|| jj=1 ||Ar|| p 
is attained at each of the 2™ points of {—1,1}" (up to scaling) and nowhere else. 
These two elements are combined in Section |4] to prove Theorcm ll.il More precisely, 
we define the matrix Z = (M T aA T ) T , where we will pick a to be a large number 
depending on n,p ensuring that the maximum of ||Zx|| p /||x|| p occurs very close to 
vectors x e { — 1,1}". As mentioned several sentences ago, the value of ||Ax|| p is 
the same for every vector x £ { — 1, 1}"; as a result, the maximum of ||Zx|| p /||x|| p is 
determined by the maximum of ||Mx|| p on {—1, 1}", which is proved in Section[2]to 
be hard to compute. We conclude with some remarks on the proof in Section [5j 

2. The || • ||oo,p norm. We now describe a simple construction which relates the 
oo,p norm to the maximum cut in a graph. 

Suppose G = ({1, . . . ,n},£7) is an undirected, connected graph. We will use 
M(G) to denote the edge- vertex incidence matrix of G; that is, M{G) € i?'^*"; we 
will think of columns of M{G) as corresponding to nodes of G and rows of M(G) as 
corresponding to the edges of G. The entries of M(G) arc as follows: orient the edges 
of G arbitrarily, and let the i'th row of M(G) have +1 in the column corresponding 
to the origin of the i'th edge, —1 in the column corresponding to the endpoint of the 
i'th edge, and at all other columns. 

Given any partition of {1, . . . , n} = S U S c , we define cut(G, S) to be the num- 
ber of edges with exactly one endpoint in S. Furthermore, we define maxcut(G) = 
maxg C { 1 n j cut(G, S). The indicator vector of a cut (S,S C ) is the vector x with 
Xi = 1 when i € S and Xi = —1 when i £ S c . We will use cut(x) for vectors 
x e { — 1,1}" to denote the value of the cut whose indicator vector is x. 

Proposition 2.1. For any p>l, 

max ||M(G)x|| p = 2maxcut(G) 1/p . 

IMU<i 



Proof. Observe that ||M(G)x|| p is a convex function of x, so that the maximum 
is achieved at the extreme points of the set ||x||oo ^ 1 5 i-e. vectors x satisfying 
Xi = ±1. Suppose we are given such a vector x; define S = {i \ Xi = 1}. Clearly, 
||M(G)x||p = 2 p cut(G, S). From this the proposition immediately follows. □ 

Next, we introduce an error term into this proposition. Define /* to be the optimal 
value /* = max|| x || (Xj <i ||M(G)x|| p ; the above proposition implies that (/*/2) p = 
maxcut(G). We want to argue that if / apP rox is close enough to /*, then (/ apP rox/2) p 
is close to maxcut(G). 
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Proposition 2.2. Ifp>l, \f* - / approx | < ef* with e < 1, then 



/ap P roxV_ maxcut(G) 



< 2 p - 1 pe ■ maxcut (G) 



Proof. By Proposition 12.11 maxcut (G) = (f*/2) p . Using the inequality 
\a p -b p \ < \a-b\pmax(\a\,\b\y-\ 

wc obtain 



fa 



maxcut (G) 



< 7; If* ~ Japproxlpmax IV'' 



r u 

2 ' 2 



p-i 



It follows from e < 1 that / apP rox < 2/*. We have therefore 

J approx \ 



— maxcut (G) 



^^Ir-Zapproxl-P-C/T-^IK/T, 



where wc have used the assumption that \ f* — /approx | < e/*- The result follows then 
from maxcut(G) = (/*/2)p. □ 

We now put together the previous two propositions to prove that approximating 
the || ■ ||oo,p norm within some fixed relative error is NP-hard. 

Theorem 2.3. For any rational p > 1, and S > 0. unless P = NP, there is 
no algorithm which given a matrix with entries in {—1,0,1} computes its p- norm to 

relative error e = ((33 + S)p2 p ' 1 ) with running time polynomial in the dimensions 
of the matrix. 



Proof. Suppose there was such an algorithm. Call /* its output on the \E\ X n 
matrix M(G) for a given connected graph G on n vertices. It follows from Proposition 
2~2lthat 



approx 



maxcut (G) 



< 



2 p - 1 p 1 

■— j-maxcut(G) = ; — jmaxcut(G). 



(33 + 8)p2 



33 + 5 



Observing that 
32 + 5 



33 — (~ S I 1 
34 + (5 maxcut ( G ) = M + 5 ( maxcut ( G ) ~ 33 + ^ maxcutCG) 



the former inequality implies 
32 + 5 



34 + 5 



maxcut (G) < 



33 + 5 / / app 

34 + 5 



< maxcut (G). 



Since p is rational, one can compute in polynomial time a lower bound V for §jr§ (/approx/ 2) p 
sufficiently accurate so that V > maxcut (G) > j|maxcut(G). However, it has 



been established in [5] that unless P = NP, for any 8' > 0, there is no algorithm 
producing a quantity V in polynomial time in n such that 



+ maxcut(G) < V < maxcut(G). 



□ 



Remark: Observe that the matrix M(G) is not square. If one desires to prove 
hardness of computing the oo,p-norm for square matrices, one can simply add \E\ — n 
zeros to every row of M(G). The resulting matrix has the same oo,p norm as M(G), 
is square, and its dimensions are at most n 2 x n 2 . 

3. A discrete set of exponential size. Let us now fix n and a rational p > 2. 
We denote by X the set { — 1, 1}™, and use S(a,r) = {x e R n | \\x — a\\ p — r} to stand 
for the sphere of radius r around a in the p-norm. We consider the following matrix 
in R 2nxn : 



A 




and show that the maximum of for x £ S(0,n 1 / p ) is attained at the 2" vectors 

in X and no other points. For this, we will need the following lemma. 

Lemma 3.1. For any real numbers x,y and p > 2 

\x + y\* + \x-y\> <2^ 1 (\x\" + \y\"). 
In fact, \x + y\ p + \x — y\ p is upper bounded by 



2 V~l W + W) _ ( p(p _ 1)| N + \y\\ p - 2 



2 Marl 



\y\ 



,p-2 



where the last term on the right is always nonnegative. 

Proof. By symmetry we can assume that x > y > 0. In that case, we need to 
prove 

(x + y) p + (x- y) p < 2 p ~ 1 (x p + y p ) - {x ~^ (p(p - l)(x + y)^ 2 - 2(x - y)^ 2 ) . 
Divide both sides by (x + y) p and change variables to z = (x — y)/(x + y): 

l + Z p < + ^ + _ (PjP^l z 2 _ V 

2 V 4 2 



The original inequality holds if this inequality holds for z £ [0, 1]. Let's simplify: 

2 + z p < (1 + zf + (1 - zf - P<yP ~ l \ 2 . 

Observe that we have equality when z — 0, so it suffices to show the right-hand side 
grows faster than the left-hand side, namely: 

z p - x < (1 + zf- 1 - (1 - zf- 1 -(p- l)z, 

and this follows from 

(1 + z)^ 1 > 1 + (p - l)z > (1 - zf" 1 + z^ 1 + (p - l)z, 
where we have used the convexity of f(a) = a 1 '- 1 . □ 

Now we are ready to prove that vectors in X optimize ||Ax||p/||a;||p, or, equiva- 
lently, optimize over the sphere 5(0, n 1 ' p ). 

Lemma 3.2. For any p > 2 ; the supremum of \\Ax\\ p over S(0,n 1 ^ p ) is achieved 
by any vector in X . 

Proof. Observe that ||^4a;||p = n2 p for any x £ X. To prove that this is the largest 
possible value, we write 

(3.1) \\Ax\\ p p = £>i - x l+1 \ p + \xi + x l+1 \ p , 

i=l 

using the convention n + 1 = 1 for the indices. Lemma |3. II implies that 

\ Xi - x i+1 \ p + \ Xi + x l+1 \ p < 2P- 1 (\xi\ p + \x t+1 \ p ) . 
Applying this inequality to each therm of (|3.1[) and using ||x||^ = n, we obtain 



\\Ax\\l < £ ^ (M P + k+iH =2^ \ Xi \* = 2"n. 

□ 



i=l i=l 



|P 

I 

decreases as we move a little bit away from the set X while staying on S(0, n 1 ^) 



Next, we refine the previous lemma by including a bound on how fast H-Axl, 



Lemma 3.3. Let p > 2,c £ (0,1/2] and suppose y £ S(0,n 1 / p ) has the property 

that 

(3.2) min \\y - x\\oo > c. 



Then, 



Proof. We proceed as before in the proof of Lemma 13.21 until the time comes to 
apply Lemma 13.11 when we include the error term which we had previously ignored: 



\\Ay\\ p p < n2*~£(lwl - |^ +1 |) 2 (p(p - 1) || Vi | + \y i+ x\\ P 2 - 2 || yi | - 



Note that on the right-hand side, we are subtracting a sum of nonnegative terms. The 
upper bound will still hold if we subtract only one of these terms; so we conclude that 
for each k, 

\\Ay\\l <n2^~- A {\y k \ - \y k+1 \f (p(p-l) \\y k \ + \y k+1 \\ p - 2 - 2 \\ Vi \ - |y l+1 ||^ 2 ) . 

By assumption, there is at least one y k with \\y k \ — l| > c. Suppose first that 
\y k \ > 1. Then we have \y k \ > 1 + c, and there must be an yj with < 1 for 
otherwise y would not be in 5(0, n 1 ^). Similarly, if \y k \ < 1, then \y k \ < 1 — c and 
there is a j for which \yj\ > 1. In both cases, this implies the existence of an index m 
with \y m \ and |y m +i| differing by at least c/n and such that at least one of \y m \ and 
|j/ m +i| is larger than or equal to 1 — c. Therefore, 

\\Ay\\l <n2V-\- 2 fp(p - 1) ||y m | + |y„ l+1 ||^ 2 - 2 ||y m | - \y m ^ 2 

Now observe that ||j/m|-|y m +i|| < |j/m|+|y m +i|, and that |y TO | + |y m +i| > (1-c) > 1/2 
because c € (0, 1/2]. These two inequalities suffice to establish that the term in square 
brackets is at least (l/2)P~ 2 (p(p - 1) - 2) > (3/2*) (p - 2), so that 

Py|| p <n2P- 3(P ~ 2) C 2 . 
N flip - 2Pn 2 

□ 

4. Proof of Theorem 11.11 We now relate the results of the last two sections 
to the problem of the p-norm. For a suitably defined matrix Z combining A and 
M(G), we want to argue that the optimizer of ||Zx|| p /||a;||p is very close to satisfying 
\xi\ = \xj\ for every i, j. 

Proposition 4.1. Let p > 2, and G a graph on n vertices. Consider the matrix 



Z 



A 



with M(G) and A as in Sections^ and\3\ respectively. If x* is the vector at which the 
optimization problem max ieS ( ll i/p) achieves its supremum then 



■ II * II <r 1 

mm \\x — x Uo < 



xex 11 " - 4Pn 6 ' 



Proof. Suppose the conclusion is false; then using Lemma [3.31 with c = l/4 p n 6 , 
we obtain 

\\ Ax *\\ Pp < n2P _ 3 <r-2) n2P _ 3(p- 2) 

II Hp- 2P4 2 Pn 14 32Pn 14 



It follows from Proposition 12.11 that 

< 2 p maxcut(G) < 2 p n 2 , 

so that 

ll^l«-ll^l B + (^)'ll^l 8 ^-^ + 2^. 

Observe that the last term in this inequality is smaller than the previous one (in 
absolute value). Indeed, for p > 2, we have that 3/32? > (2/64) p ,p-2 > [(p-2)/p] p 
and 1/n 14 > n 2 /n 8p . We therefore have ||Zx*||^ < 2 p n. By contrast, let x be any 
vector in {-1, 1}". Then x € 5(0, n 1/p ) and 

\\Zx\\l > \\Ax\\l > 2 p n, 

which contradicts the optimality of x*. □ 

Next, we seek to translate the fact that the optimizer x* is close to X to the fact 
that the objective value ||Zx||p/||x|| p is close to the largest objective value at X. 

Proposition 4.2. Letp > 2, G a graph on n vertices, and 



.4 



Z = P-2 

V M{G) J 

If x* is the vector at which the optimization problem 

max llZxIlp 

x€S(0,nVp) 

achieves its supremum and x r is the rounded version of x* in which every component 
is rounded to the closest of —1 and 1, then 



\Zx*\\r-\\Z Xl \\* 



1 



Proof. Observe that x* is the same as the extremizer of the corresponding problem 
with Z instead of Z, so that x satisfies the conclusion of Proposition ETTl Consequently 
every component of x* is closer to one of ±1 than to the other, and so x T is well defined. 
We have: 

\\Zx*\\v - \\Zx t \\* = (64^n*y (\\Ax% - \\AxrW*) + (\\Mx*\\ p - \\Mx Y \\ p ). 

This entire quantity is nonncgative since x* is the maximum of \\Zx\\ on 5(0, n 1 ^). 
Moreover, ||Ae*||p — ||Ax r |||j is nonpositive, since by Proposition ^. 21 1 |Ae|| p achieves 
its maximum over 5(0, n x / p ) on all the elements of X . Consequently, 

\\Zx*\\* -\\Zx t \\l<\\Mx*\\l-\\Mx t \\l 
(4.1) < (\\Mx*\\ p - ||Ma; r ||p)i?max(||Ma;*||p, WMx^pf- 1 . 
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We now bound all the terms in the last equation. First, 

(4.2) ||Mx*|| p - ||Mx r || p < ||M|| 2 ||x* - x r || 2 < \\M\\ F Vn\\x* - x r 



4P n 6 ■ 



where we have used ||M(G)||i? = y/2 \E\ < n and Proposition 14 . 1 1 for the last inequal- 
ity. Now that we have a bound on the first term in Eq. (|4.ip . we proceed to the last 
term. It follows from the definition of M that 



\MxA\l < TP ■ [ "J | < 2 p n 2 



Next, we bound ||Ma;*||^. Observe that a particular case of Eq. (|4.2|) is: 



(4.3) 



\Mx* 



< \ \Mx T \ 



Moreover, observe that ||Mx r || p > 1 (the only way this does not hold is if every entry 
of x T is the same, i.e. ||Mx r || p = 0; but then Eq. (|4.3|) implies that ||Mx*|| p < 1, 
which is impossible since G has at least one edge), and so Eq. (|4.3[) implies ||Mx*|| p < 
2||Mx r || p , and so: 

||Afx*||J < 4 p n 2 . 

Thusmax(||Afx*||p, ||Mx r ||p) p < 4 p n 2 and thcreforet max(||Mx*|| p , ||Mx r || p )P- 1 < 4 p n 2 . 
Indeed, this bound is trivially valid if max(||Mx*||p, ||Mx r ||p) p < 1, and follows from 
a p_1 < a p for a > 1 otherwise. Using this bound and the inequality (|4.2[) . we finally 
obtain 



J, 1.5 

\Zx*\\ p -\\Zx r \\ p < 

I Hp II r||p - 4 p n 6 



p . 4P n 2 < 1 



Finally, let us bring it all together by arguing that if we can approximately com- 
pute the p-norm of Z , we can approximately compute the maximum cut. 

Proposition 4.3. Let p > 2. Consider a graph G on n > 2 vertices and the 
matrix 

p-2 



M(G) 



and let f* = \\Z\\ P . If 

(4.4) 

then 



|/a 



ri < 



(p - 2f 



— f p - n 

2p J approx 



132V" 8p+3 p' 
/64pn 8 ^ p ^ 



maxcut(G) 



< 



8 



Proof. Observe that n,p f* = max ieS ( n 1 /*) ll-^ x llp- It follows thus from Proposi- 
tion [O that 



(4.5) 



nf* p — max ||Za;||£ 



< 



1 



Recall that \\Zx\\p = \\Mx\\p + (64^n 8 ) P \\Ax\\p, and that \\Ax\\p = n2P for evei T 
Therefore, 



max||Zx||? = 



/64prf \ ^ + max || Ma .|.p = f6±P]f \ n2P 2 Pmaxcut(G), 
\p — 2 J xex " p \p — 2 J 



\p-2 

and combining the last two equations we have 

/64pn s ^ p 



(4.6) 



71 

— f* p - n 
2p 



P 



— maxcut(G) 



< 



1 



2Pn 2 



Let us now evaluate the error introduced by the approximation / apP rox 
/ 64pn 8 x p 



2^ J approx 



p-2 



maxcut(G) 



<_L + H|p _ f*P\ 
— 2Pn 2 2p > Ja -PP IOK J I 

~ 2Pn 2 ^ 2P ^ approx _ ^ |P max (/*i /approx) P 1 - 



It remains to bound the last term of this inequality. First, we use the fact that /* > 1 
and Eq. (|4.6p to argue 



(4.7) f*<-P-V < f*P < 2 P 



( 64pn s 



8\ P 



2» 



1 



— maxcut(G) + — < 2 P 



8\P 



\p-2 



where we have used maxcut(G) < n 2 and 1 < p/(p — 2) for the last inequality. By 
assumption, |/ apP rox — f*\ < 1 and since /* > 1, 



/^pr 1 o ) x <(2rr 1 <(2r) p <4P (® 



66pn J 



8\ P 



- 2 



Putting it all together and using the bound on | /approx — f*\, we obtain (assuming 
n > 1) 



— f p -n 

Qp J approx 



\~P~2 



8\P 



— maxcut(G) 



< 



1 



(p-2)P 



-2 p np 



66pn 



8\P 



< 



< 



2Pn 2 l32 p pP7i s P+ 3 p \p-2 

1 1 
2Prf + ^2 
1 



Proposition 4.4. Fix a rational p <G [l,oo) with p ^ 1,2. Unless P = NP, 
there is no algorithm which, given input e > and a matrix Z, computes \ \Z\\ P to a 
relative accuracy e, in time which is polynomial in 1/e, the dimensions of Z , and the 
bit-size of the entries of Z . 
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Proof. Suppose first that p > 2. We show that such an algorithm could be used 
to build a polynomial-time algorithm solving the maximum cut problem. For a graph 
G on n vertices, fix 



l 32 f ( \ n 8p+3 p I • ( 132 



build the matrix Z as in Proposition ^. 31 and compute the norm of Z; let / apP rox be 
the output of the algorithm. Observe that By Eq. (|4.7[) 

lizil. < 132! " 8 



so. 



/approx 1 1 Z 1 1 p 



p-2 



% < e ( 132^n 8 j < ^132" ) n^p 



It follows then from Proposition 14. 31 that 

/ approx V n / g4 ( P \ J 



p-2 

is an approximation of the maximum cut with an additive error at most 1/n. Once 
we have /approx, we can approximate this number in polynomial time to an additive 
accuracy of 1/4. This gives an additive error 1/4 + 1/n approximation algorithm for 
maximum cut, and since the maximum cut is always an integer, this means we can 
compute it exactly when n > 4. However, maximum cut is an NP-hard problem pQ. 

For the case of p <G (1,2), NP-hardness follows from the analysis of the case of 
p > 2 since for any matrix Z, \\Z\\ P = \\Z T \\ p / where l/p + 1/p' = 1. □ 



Remark: In contrast to Theorem l2.3l which proves the NP-hardness of computing the 
matrix oo, k norm to relative accuracy e = 1/C(p), for some function C(p), Proposi- 
tion H3] proves the NP-hardness of computing the p-norm to accuracy 1 / C (p)n Sp+11 , 
for some function C'{p). In the latter case, e depends on n. 



Our final theorem demonstrates that the p-norm is still hard to compute when 
restricted to matrices with entries in {—1,0, 1}. 



Theorem 4.5. Fix a rational p € [1, oo) with p ^ 1, 2. Unless P = NP, there is 
no algorithm which, given input e and a matrix M with entries in {—1, 0, 1}, computes 
||Af|| p to relative accuracy e, in time which is polynomial in e _1 and the dimensions 
of the matrix. 



Proof. As before, it suffices to prove the theorem for the case of p > 2; the case 
of p e (1,2) follows because H^Hp = H-Z^Hp' where l/p+ 1/p' = 1. 
Define 




where [•] refers to rounding up to the closest integer. Observe that by an argument 
similar to the proof of the previous proposition, computing ||Z*|| P to an accuracy 
e = (C(p)n 8p+11 )~ 1 is NP-hard for some function C(p). But if we define 

/ A \ 

A 



where A is repeated (^^^2 n8 ^j 



A 

V M J 



times, then 



\\Z**\\ P = ||Z*|| P . 

The matrix Z** has entries in { — 1,0, 1} and its size is polynomial in n, so it follows 
that it is NP-hard to compute ||i?**|| p within the same e. □ 

Remark: Observe that the argument also suffices to show that computing the p-norm 
of square matrices with entries in { — 1, 0, 1} is NP-hard: simply pad each row of Z** 
with enough zeros to make it square. Note that this trick was also used in Section [2] 

5. Concluding remarks. We have proved the NP-hardness of computing the 



matrix p-norm approximately with relative error e = 1/C(p): 



,8p+ll 



, where C(p) is 



some function of p; and the NP-hardness of computing the matrix oo,p norm to some 
fixed relative accuracy depending on p. We finish with some technical remarks about 
various possible extensions of the theorem: 



Due to the linear property of the norm ||<xA|| = \a\ \\A\\, our results also 
imply the NP-hardness of approximating the matrix p-norm with any fixed 
or polynomially growing additive error. 

Our construction is also implies the hardness of computing the matrix p-norm 
for any irrational number p > 1 for which a polynomial time algorithm to 
approximate x p is available. 

Our construction may also be used to provide a new proof of the NP-hardness 
of the || • ||p j9 norm when p>q, which has been established in [5]. Indeed, it 
rests on the matrix A with the property that max ||Ax|| p /||x|| p occurs at 
the vectors x E { — 1,1}™. We use this matrix A to construct the matrix 
Z = {aA M) t for large a, and argue that max||Zx|| p /||x|| p occurs close to 
the vectors x € { — 1, 1}™- At these vectors, it happens Ax is a constant, so 
we are effectively maximizing ||A/x|| p , which is hard as shown in Section [SJ 
If one could come up with such a matrix for the case of the mixed 1 1 ■ | | Pi9 norm, 
one could prove NP-hardness by following the same argument. However, when 
p > q, actually the very same matrix A works. Indeed, one could simply argue 
that 



\Ax\ 



- max ■ 



max 



\Ax\\ q \\x\ 



and since the maximum of ||a;|| g /||x||p when 1 < q < p < oo occurs at the 
vectors x G { — 1, 1}™, we have that both terms on the right are maximized at 
x =E { — 1,1}", that is where ||j4x|| g /||x|| p is maximized. 
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• Finally, we note that our goal was only to show existence of a polynomial-time 
reduction from the maximum cut problem to the problem of matrix p-norm 
computation. It is possible that more economical reductions which scale more 
gracefully with n and p exist. 
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